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ABSTRACT 

We calculate the 0(arrq /rriyy) corrections arising from diagrams involving the top- 
quark loops to the light neutral Higgs boson production via qq' — > WH at the 
Fermilab Tevatron in both the standard model and the minimal supersymmetric 
model. In contrast to the QCD correction which increases the tree-level cross sec- 
tion, the corrections imply a few percent reduction in the production cross section 
relative to the tree-leve results. 



PACS number: 14.80.Bn, 14.80. Cp, 13.85. QK, 12.60. Jv 



I. INTRODUCTION 



The Higgs boson is the only particle of the Standard Model(SM) which has not been 
discovered so far. The direct search in the LEP experiments via the e + e~ — > Z*H yields a 
lower bound of ~ 77. 1 GeV on the Higgs mass Jl| . This search is being extended at present 
LEP2 experiments, which will explore up to a Higgs boson mass of about 95 GeV via 
e + e~ — > ZH by the year 2000 0. Much higher masses will be explored by the CERN Large 
Hadron Collider (LHC). Since it will be some years before the LHC comes into operation it 
is worth considering wherther the Higgs boson can be discovered from the existing hadron 
collider, the Tevatron. Much study has been made in the detection of a Higgs boson at 
the Tevatron . It was recently pointed out || that a light Higgs boson of mass 60 GeV 
< rriH < 130 GeV can be observable at the Tevatron with CM energy a/s = 2 TeV and 
sufficient integrated luminosity, 30— 100/6 -1 , through the production subprocess qq' —>■ WH, 
followed by W — > iv and H — > bb. Since the expected number of events is small, it is 
important to calculate the cross section as accurately as possible. In Ref. the 0(a s ) 
QCD corrction to the total cross section to this process have been calculated, and the QCD 
correction were found to be about 12% in the MS scheme at the Fermilab Tevatron and 
the LHC in the SM. In general, the SM electroweak corrections are small comparing with 
the QCD correction. Beyond the SM, the electroweak corrections might be enhanced, since 
more Higgs bosons with stronger couplings to top or bottom quarks are involved in some 
new physics models; for example, the minimal supersymmtric model(MSSM) 0] J7|, which 
predict that the lightest Higgs boson ho be less than 140GeV. Therefore, it is worthwhile 
to calculate the electroweak corrections to the light Higgs boson production via qq' — > Wh Q . 
In this paper we present the calculation of the top quark loop correction of order avr? t / m^y 
to the Higgs boson production at the Fermilab Tevatron in both the SM and the MSSM. 
These corrections arise from the virtual effects of the third family (top and bottom) of quark, 
neutral and charged Higgs bosons. And we shall present the complete calculations of the 
electroweak radiative corrections to this process in a future publication || . 
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II. CACULATIONS 



The Feynman diagrams for the lightest Higgs boson production via q(pi)q'{p 2 ) — > 
W{ki)ho{k 2 ) , which include the top quark loop corrections of order amf/m^y to the pro- 
cess qq' — > Who, are shown in Fig 1. We perform the calculation in the 't Hooft-Feynman 
gauge and use dimensional regularization to all the ultraviolet divergences in the virtual 
loop corrections utilizing the on-mass-shell renormalization 0, in which the fine-structure 
constant a and the physical masses are chosen to be the renormalized paramenters, and the 
finite parts of the countertems are fixed by the renormalization conditions. As far as the 
parameters j3 and a, for the MSSM we are considering, they have to be renormalized, too. 
In the MSSM they are not independent. Neverthless, we follow the approach of Mendez 



and Pomarol [10| in which they consider them as independent renormalized parameters and 



fixed the correspoding renormalization constant by a renormalization condition that the 



on- mass-shell H + £v£ and h§ll couplings keep the forms of Eq.(3) of Ref. [10 to all order of 
perturbation theory. 

We define the Mandelstam variables as 

s = (P1+P2) 2 = (h + k 2 ) 2 
i = (Pi ~ h? = (P2 - k 2 f 

u = (pi- k 2 f = (p 2 - h) 2 . (1) 
The relevant renormalization constants are defined as 

m wo = m w + 5m wi m zo =rn 2 z + 5m 2 z , (2) 

tan/3 = (1 + 5Zp) tan/3, sina = (1 + SZ a ) sin a, (3) 

Wf» = Zfw±» + iZ^PH*, Ht = (1 + SZ^H^, (4) 

ho = (l + 5h Q y' 2 h + Z)£ H H, H = (1 + 5Z H ) l ' 2 H + Z^h (5) 

Taking into account the 0(am 2 /m 2 v ) corrections, the renormalized amplitude for qq' — > 
Who can be witten as 



M ren = M + 8M self + SM 



vertex 



(6) 



where M is the amplitude at the tree level, 8M sel f and SM vertex represent the corrections 
arising from the self-energy and vertex diagrams, respectively. M n is given by 

e 2 m w sin(/3 — a 



V2{ 



rn 



w 



s) sin 9 2 



d(p 2 )<fP L u(p 



(7) 



where Pl,r = (1 T 7s)/2- 5M sel f is given by 
6M sei f = 5m 2 w + (m 2 w - s)5Z w 



s - mfc 

-A 



+ 



N c e mw sin(/3 — at) 



6sm 2 - 2s 2 + 3m 2 {m\ - 2s) 



288 v / 2V 2 s(-m 2 v + s) 2 sin 9* 
x B (0, ml, m 2 t ) + 3(— m\ - sm 2 t + 2s 2 )B (s, ml, m\) d(p 2 )iPLu{pi) 



(8) 



with 



5m 2 w 



N c e 2 m 2 
967T 2 sin # 2 



-2 + 2B (0, m 2 b , m 2 ) - B (m 2 w , m 2 b , m 2 ) - 45 (0, m 2 , m 2 ] 



+ 



m. 



in 



w 



B (0, ml, m 2 ) - B {m 2 w , ml, m 2 ) 



(9) 



5Z 



N r e 2 



w 



288m^7r 2 sin 9 2 w 



2m w + 3m t B o (0, ml, m\) - 3{m\ + 2m\ v ) 



x B (m w ,m b ,m t ) + 3m w (m t + m t m w - 2m w )G(m w ,m b ,m 



(10) 



Here and below, Bq, Cq, Cj and Cjj is the two-point and three-point scalar integrals, defi- 



nitions for which can be found in Ref. |TT[ and G is the derivative of Bo which is expressed 
as 



G(M 2 ,Mf,M 2 
5 M vertex is given by 

5NI ver ^ ex 



2 n#2 , dB (k 2 ,M 2 ,M 2 ) 



dk 2 



\k 2 =M 2 - 



(11) 



1 jmgg - 8m% 

2 dZh0 + 2(m 2 z -m 2 w ) 



5m^ 5m 2 w 

W 



in 



+ cot(/? — a)(Z]H + sin /3 cos /?5Zg — tana5Z a 

+ fr tex d(P2W L u( Pl ) 

+ fr teX AV2)hPLu{p x )e.p x 
+ fr^dip^PMpJe.pi, 



(12) 
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N c e 2 m 2 
3 cos 61,-n 2 



B (0,m 2 ,m 2 ) B (0,m 2 t ,m 2 ) 2sm6 2 w B (0,m 2 ,m 2 ) 



6 



16 sin 9 2 
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B (m 2 z ,m 2 ,m 2 ) B (m 2 z ,m 2 ,m 2 t ) 2sm9lB ( y m 2 z ,m 2 ,m 



32 sin 6 2 w 



+ 



8Z ho = 



-N c e 2 cos 2 (a) csc 2 ((3)m 2 
32m 2 y 7r 2 sin 6^ 



B (m ho ,m t ,m t ) 



+ (m ho -4m t )G(m ho ,m t ,m t ) 



7 l/2 
J Hh 



N c e m 2 cos(a) esc ((3) sin(o;) 
32 [mil 2 — nrih, 2 ) mw 2 n 2 sin6*„ 
+ (mL - Am 2 )B (m 2 h m\, m 2 ) 



-2m 2 t -2m 2 t B (0,m 2 t ,m 2 ) 



6Z* 



5m 2 7 — 5m 2 v 5m 2 ? bm 2 M ' 



+ 



w 



2(to| - m 2 w ) 2m\ 2m 2 w 2 



K ry m W ryl/2 

-0Zj H ± — -Zj wh± , 



tan j3 



sin 2 pSZ fi + 5m \- 5m k - H + ^ - ^Z,. — -Z 



cos a 



2(m| - m 2 w ) 2m| 2to^ 2 



sin a 



1/2 
fZ7i > 



8Z H ± 



N c g 2 m 2 
32m 2 v ir 2 



- B (m 2 H± , m 2 b , to 2 ) cot 2 ((3) 



+ (m 2 - m 2 H ±) cot 2 (/3)G(m 2 H ±, m 2 , m 2 



1/2 N c g 2 m 2 t cot(p) r 
H±iy 32to 2 , ± to^tt 2 



m 2 B (0,m 2 b ,m 



+ ( m #± -m 2 )B (m 2 H ±,ml,m 2 



n 



vertex 



-N c e A m 2 



32 \/2 tovk 7r 2 (mvi/ 2 — s) sin 0„ 



-2B (s,ml,m 2 t ) 



+ (- 2m t ~ t) C o( m l i m wi s, m\, m\, m 2 b ) 
+ {-2m 2 w - i)Ci(m 2 w , s, ml, m 2 t , m 2 b , m\) 



+ (- m l ~ ^i)C 2 {m 2 w , s, m 2 ho ,m 2 t , m 2 b , m 2 ) 
+ AC m (m 2 w ,s, m 2 ho ,m 2 t ,m 2 b ,m 2 )] , 



fr tex = ~77~~/k 2 N ( C ^"f . fl 4 [-C {ml,m^mlmlml) 

IbyZmw^ [rriw s) smu w L 

n /_2 2 2 2 2\ 

- Ci(m w , s,m ho ,m t ,m b ,m t ) 

- 3C 2 (m^, §, m 2 hQ , m 2 , m 2 b , m\) 

nr^i ! 2 - 2 2 2 2\ 

- 2C l2 {m w ,s,m ho ,m t ,m h ,m t ) 

- 2C 22 (m 2 w , s, m 2 ho} m 2 } to 2 , to 2 )] , (21) 

pveri.PT, N c 6 771+ r . n o 2 2 2\ 

/ 3 = TT"^ , ^ . 7 4 ^K, s, r<, to; ?, to 2 , to 2 ) 

o/i / 2 - 2 2 2 2\ 

- 2C l2 {m w ,s,m hQ ,m t ,m b ,m t ) 

- 2C 22 (m 2 w , §, m 2 ho ,m 2 , m 2 b , m 2 )\ . (22) 
The corresponding amplitude squared for the process qq' — > Wh can be written as 

J2 \M ren \ 2 = £ |M | 2 + 2ReJ2{5M self + bM vertex )M\, (23) 

where the bar over the summation recalls average over initial partons spins. The cross 
section of process qq' — > Wh is 



with 



a 



u = < + < - ; _ V(j _ (mfci + mw)2)(j _ (mto _ mw)2)/2 

imax = — h V ( 5 ~ ( mft o + m vy) 2 )(s - {m ho - m w ) 2 )/2. (25) 

The total cross section of PP — > qq 1 — >■ VF/io can be obtained by folding the a with the 
parton luminosity 

a(s) = f 1 dz^a(qq -> W/i at s = ^ 2 s), (26) 

where a/s and a/I is the CM energy of PP and g<f, respectively, and dL/dz is the paron 
luminosity, which is defined as 

g = 2zJ^f q/P (x,q 2 )f g//P (z 2 /x,q 2 ), (27) 

where f q /p(x,q 2 ) and f q , /p(z 2 /x,q 2 ) are the parton distribution function fll2"| . 



167TS 2 

spins 



E 1^1 ^ (24) 
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III. NUMERICAL RESULTS 



In the following we present some numerical results. In our numerical calculations, the 
SM parameters were taken to be mw = S0.33GeV, mz = 91.187GeV, m t = 176GeV, 
nib = 4.5GeV and a(mw) = j^g- Moreover, we use the relation |7j between the Higgs 
boson masses m ho H A H ± and parameters a, (3 at one-loop, and choose rrih and tan/3 as two 



independent input parameters. As explained in Ref. [10], there is a small inconsistency in 
doing so since the parameters a and (3 of Ref. [7[ are not the ones defined by the conditions 



given by Eq.(3) of Ref. |T0[. Nevertheless, this difference would only induce a higher order 
change |10| . We will limit the value of tan/5 to be in the ranges 2 < tan/3 < 30, which are 
consistent with ones required by the most popular MSSM model with scenarios motivated 
by current low energy data (including a s , Rb and the branching ratio of b — > 37). 

In Fig. 2 we present the tree-level total cross sections versus the Higgs boson mass in 
both the SM and the MSSM for the different values of tan/3, using the CTEQ3L parton 



distributions ||12|| . Figure 2 shows the total cross sections in the SM always are larger than 
ones in the MSSM, and they are almost same only for tan/3 = 2, or the mass of the Higgs 
boson approch to 130 GeV. 

In Fig. 3 we shows the top quark loop corrections of order amf/m^ to the total cross 
sections. From Fig. 3 one sees that in the SM the corrections are not sensitive to the mass of 
the Higgs boson and amounts to 1% ~ 2% reduction in the cross section. And in the MSSM 
the corrections depend strongly on the values of rrih for all tan (3. Especially for tan (3 > 2, 
such correction can reach about —4% when = 60GeV, but the correction is only about 
— 1% if rrih = !30GeV. Since QCD correction increases the tree- level total cross sections 
by about 12% it is necessary for an accurate calculation of the cross sections to include 
the top quark loop corrections, which typically imply a few percent reduction in the cross 
sections. 

In conclusion, we have calculated that the top quark loop corrections of order avr^jm^ 
to the neutral Higgs boson production via qq' — > WH at the Fermilab Tevatron in the SM 
and the MSSM. In contrast to the QCD corrections, such corrections reduce the tree-level 
total cross sections by about 1% ~ 2% in the SM, and 1% ~ 4% in the MSSM. 
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FIG. 2. Tree- level cross sections as a function of the Higgs boson mass of the process qq' — ► Who 
with yO? = 2TeV at Tevatron. 
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FIG. 3. Relative corrections Sa/ao as a function of the Higgs boson mass of the process 
qq' — > Who with y/s = 2TeV at Tevatron. 
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